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ABSTRACT 

We have analyzed H and if^-band images of the Arches cluster obtained using the NIRC2 instru- 
ment on Keck with the laser guide star adaptive optics (LGS AO) system. With the help of the LGS 
AO system, we were able to obtain the deepest ever photometry for this cluster and its neighborhood, 
and derive the background-subtracted present-day mass function (PDMF) down to 1.3 M0 for the 
5"-9" annulus of the cluster. We find that the previously reported turnover at 6 Mq is simply due 
to a local bump in the mass function (MF), and that the MF continues to increase down to our 
50 % completeness limit (1.3 M©) with a power-law exponent of F = —0.91 for the mass range of 
1.3 < A// Mq < 50. Our numerical calculations for the evolution of the Arches cluster show that the 
F values for our annulus increase by 0.1-0.2 during the lifetime of the cluster, and thus suggest that 
the Arches cluster initially had F of —1.0 ^ —1.1, which is only slightly shallower than the Salpeter 
value. 

Subject headings: Galaxy: center — open clusters and stellar associations: general — stars: luminosity 
function, mass function — celestial mechanics, stellar dynamics 



1. INTRODUCTION 

The stellar initial mass function (IMF) is the most pri- 
mary product of the star formation process. With this 
one relation, it is possible to directly probe the crucial 
predictions of star formation models. Surprisingly, the 
IMF is approximately universal, following the Salpeter 
law (Salpeter 1995) for masses between 1 and 120 M© 
(Kroupa 2002). On the low mass end, there appears to 
be a universal rollover around 0.8 M©, and the high mass 
end seems to be truncated by a sharp cutoff near 150 M© 
(Weidner & Kroupa 2004; Figer 2005). 

However, star formation theories predict that the lower 
mass cutoff (mi) should be a function of environmental 
parameters, i.e. magnetic field strength and cloud tem- 
perature (e.g., Bonnell, Larson, & Zinnecker 2006). In 
the extreme environment of the Galactic center (GC), 
some of these models predict an elevated lower mass cut- 
off/rollover with respect to that observed in the disk 
(Morris 1993). If this prediction is true, then the GC 
should have an abnormally bright lower mass cutoff, 
something that has not been directly observed anywhere 
in the Universe. 

The best places to trace the IMF near the GC are the 
two young star clusters therein: the Arches and Quintu- 
plet clusters (see Figer et al. 1999 for findings and early 
studies on these clusters) . These clusters are very young 
(2-4 Myr), compact (< 1 pc), and only 20-30 parsecs 
away from the GC in projection, while they appear to 
be as massive as the smallest Galactic globular clusters 
(- 2 X 10* M©; Kim et al. 2000). Of the two clusters, 
the Arches is preferred for studies of the low-end mass 
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function (MF) as the Quintuplet is significantly more 
dispersed. 

While the Arches cluster's compactness makes its stars 
easily identifiable for their high spatial density projected 
onto the background population, this same property 
makes confusion a problem at faint magnitudes. Note 
for scale that a 1 M© dwarf has K ~ 21 in the cluster. 
For instance, the luminosity function (LF) obtained from 
Hubble Space Telescope (HST) NICMOS observations by 
Figer et al. (1999) for the whole cluster was 50 % com- 
plete at 6M0. The adaptive optics (AO) science demon- 
stration data observed with the Gemini/Hokupa'a sys- 
tem were analyzed by Yang et al. (2002) and Stolte et al. 
(2002), who obtained a 50 % completeness level at 5.6 M© 
and a 75 % completeness level at 6 M© for the whole clus- 
ter, respectively. Stolte et al. (2005) have analyzed the 
natural guide star (NGS) AO observations of the cluster 
with the Very Large Telescope NAOS/CONICA system, 
and derived the MF of the intermediate region of the 
cluster (5"-10") down to 2.3 M© at a 50 % completeness 
limit. 

Although the last study reached the faintest magni- 
tude to date, due to the lack of suitable NGS candi- 
dates for neighboring fields, it lacks the observations of 
nearby background populations, which should be sub- 
tracted from the LF/MF of the cluster field to obtain 
the true cluster LF/MF. It had no choice but to com- 
pare the MF from the central and intermediate regions 
of the cluster (r < 10") with that from the outer region 
(r > 13.4") to infer the MF of the cluster only. However, 
while there are almost no bright cluster stars in this outer 
region, a non-negligible number of faint cluster stars may 
be present there due to mass segregation. For example, 
the comparison between HST observations and TV-body 
simulations by Kim et al. (2000) indicates that the tidal 
radius of the Arches cluster is 25"-30". Thus a reliable 
estimation of the LF/MF of the nearby background pop- 
ulations would require observations of control fields that 
are at least 30" away from the cluster center. 

For this reason, we have carried out Keck/NIRC2 laser 
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Fig. 1. — Locations of our images in the galactic coordinates. 
The control fields, which were taken to estimate the stellar density 
of the nearby background population, are separated ~ 60" from 
the cluster center. 

guide star (LGS) AO observations of the Arches cluster 
and its nearby background populations to take advantage 
of the LGS AO system that does not have the distance 
limit from NGS candidates. Here, we present the analy- 
sis of these observations along with results from several 
relevant A'^-body simulations. 

2. OBSERVATIONS AND DATA REDUCTION 

The images were obtained using the Keck/NIRC2 LGS 
AO system on 2006 May 4 & 5. The Arches cluster 
(a = 17'^45'"50".59, 5 = -28°49'20".3; J2000) was ob- 
served with the medium (20" x 20", 0".0198 pixel"!) and 
wide (40"x40", 0".0397 pixel"!) field cameras, and three 
nearby control fields, separated ^ 60" from the cluster 
center, were imaged with the medium field camera in or- 
der to sample the background population. The locations 
of these fields in galactic coordinates are shown in Fig- 
ure [TJ All fields were imaged in H and Kg bands in a 
9-point dither pattern with a leg size of 2". The multiple 
correlated double sampling (MCDS) mode was used with 
10 coadds and a 10 s integration time per coadd, giving 
a total exposure time of 900 s per field. Obtained Strehl 
ratios for the {H) images are 0.17 (0.08), 0.19 (0.08), 
0.25 (0.06), 0.13 (0.06), and 0.12 (0.06) for the medium 
cluster field, control fields A, B, C, and wide cluster field, 
respectively, resulting in FWHMs of 70-100 mas for Ks, 
and 61-98 mas for H . 

For data reduction, the dithered images were first com- 
bined by maximizing the cross-correlation. Star finding, 
point-spread function (PSF) building, and PSF fitting 
procedures were performed using the DAOPHOT pack- 
age within the Image Reduction and Analysis Facility 
(IRAF). The medium field camera PSFs were moder- 
ately elongated toward the center of the image due to a 
relatively small isoplanatic patch obtained during the ob- 
servations, thus we have allowed the PSFs to be quadrat- 
ically variable along both axes of the field. 

Typical DAOPHOT magnitude errors were less than 
0.05 mag down to Kg — 19 for stars in the cluster field 
with r > 5" {Kg — 18 for the wide field camera) and 
Ks ~ 19.5 for stars in the nearby fields. We do not 
include stars with errors greater than 0.2 mag {Kg ^ 
21.5) in our analyses. The elongation of the PSF of the 
medium camera was significant at the four corners of the 
images, so we excluded stars farther than 10" from the 
image center. 

Photometric calibration was done against the 
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Fig. 2. — Dereddened color-magnitude diagram of the stars with 
DAOPHOT magnitude errors less than 0.2 mag in the 5"-9" an- 
nulus of the cluster field (left) and in control fields A, B, and C. 
Also shown are the solar-metallicity isochrones for 2 Myr ( left) and 
1 Gyr {right) from the Geneva set of models. The limiting mag- 
nitude in J? is ~ 1 mag brighter than Ks , so the number of stars 
shown in the left panel is much smaller than that used for deriving 
our mass function later, which is based on the Xs-band photometry 
only. 

HST/mCMOS photometry by Figer et al. (1999) 
assuming that H and Kg magnitudes are identical to 
NICMOS F160W and F205W magnitudes for the stars 
in our fields (Kim et al. 2005 and Stolte et al. 2002 
indicate that the transformation from H and K' to 
F160W and F2Q5W filters has a very weak dependence 
on color for the stars in the Arches cluster). 

Completeness tests were performed for all images by 
adding artificial stars to the observed images with the 
PSF obtained from our reduction. 800 artificial stars 
were used per magnitude bin (a total of 7200 stars), 
but only 200 stars were added per test image. For our 
analyses in the following sections, we will only consider 
magnitude and mass bins with a completeness fraction 
larger than 50 %. Due to the relatively irregular shape 
of the AO PSFs, the DAOPHOT star-finding algorithm 
resulted in many bogus stars around bright stars, and we 
have excluded from our photometry the faint stars that 
are ~ 4 mag fainter than a nearby (within ~ 0.55" for the 
medium field camera and ~ 0.8" for the wide field cam- 
era), brighter star. Since the central region of the cluster 
is dominated by bright stars, we will analyze the LF/MF 
of the intermediate region (5"-9") where the photometry 
of the faint stars is limited by the background confusion. 
But there are still a few bright stars in the intermediate 
region that can hamper the recovery of nearby fainter 
stars. Since we are more interested in the incomplete- 
ness of faint stars by background confusion than that by 
nearby bright stars, we exclude artificial stars that are 
^ 4 mag fainter and within ~ 0.55" 0.8" for the wide 
field) from a nearby star from calculating the complete- 
ness fractions. 

3. RESULTS 

3.1. Extinction 

For our analyses, we adopt the extinction law by Rieke, 
Ricke, & Paul (1989; Ah = 1.56Ak) and the distance 
modulus for the GC of 14.52 mag (8 kpc). Stolte et al. 
(2002) reported color variation of stars as a function of 
the distance from the cluster center and interpreted it as 
a local depletion of dust (and thus a smaller extinction) 
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Fig. 3. — Dereddened, completeness corrected Ks-band lumi- 
nosity function for control fields A, B, and C, and the outer region 
(r > 15") of the cluster field from the wide field camera, scaled 
to the area of 5"-9" annulus. The histograms are shown only for 
magnitude bins whose recovery fraction is more than 50 %. Also 
shown are the Poisson error bars for the control field A. 

at the cluster center due to winds from massive stars 
or photo-evaporation of dust grains by the intense UV 
radiation field. Our data confirm such variation, and we 
determine the extinction values for the 5"-9" annulus 
of the cluster field, the outer (r > 15") region of the 
cluster field, and three control fields separately. Then we 
assume that stars in the same field or annulus have the 
same extinction value. 

We compare the 2 Myr, solar-metallicity isochrone 
from the Geneva set of models (Lejeune & Schaerer 2001) 
to our color-magnitude [H — Kg vs. Ks) diagram (Fig. [2]) 
and obtain an average extinction at Kg of 3.1 ± 0.19 mag 
for the 5"-9" annulus of the cluster field, which is con- 
sistent with the extinctin value estimated by Stolte et al. 
(2005) for the same region. 

Our control fields and outer region of the cluster field 
are expected to be composed of stars with a variety of 
ages and metallicities (Figer et al. 2004). However, 
H — Ks colors of the stars with intrinsic Kg magni- 
tudes of 1-3 mag, which have undereddened, apparent 
Ks magnitudes of 18-20 mag at the GC, vary only lit- 
tle (< 0.05 mag; see Kim et al. 2005, for example) 
for various ages and metallicities. So we compare the 
1 Gyr, solar-metallicity isochrone to the H — Ks color 
distribution of stars with underredened, apparent Ks of 
18-20 mag, and find that our control fields A, B, & C 
have Ak, of 2.9 ± 0.22, 3.1 ± 0.17, & 2.5 ± 0.25 mag, 
respectively, and the outer region of the cluster field has 
Ak^ = 3.0 ± 0.25 mag (quoted uncertainties are the av- 
erages of the extinction spread). Since the amounts of 
extinction spread are similar in our cluster and control 
fields, the differential extinction in each field will not 
cause a problem to our subtraction of background popu- 
lation from the cluster field. 

3.2. Background Population 

Figure [3] compares the dereddened, completeness- 
corrected i^g-band LFs for the control fields and that 
for the outer region (r > 15") of the cluster field ob- 
tained with the wide field camera. While the LFs of con- 
trol fields A and C are quite similar, control field B and 
the outer region of the cluster field generally have higher 
LFs than the former. This is probably because the con- 



FlG. 4. — Dereddened Xs-band luminosity functions for the 
5"-9" annulus of the cluster field (solid) and the average of the 
control fields A, B, and C (dotted) scaled for the area of 5"-9" 
annulus. Thin lines are for the raw LFs and thick lines are for the 
completeness-corrected LFs. The Poisson error bars are shown for 
the latter. The arrow indicates the 50 % completeness limit for the 
5"-9" annulus of the cluster field. 

trol field B is located closer to the Galactic plane where 
the stellar density is larger and a part of the outskirt and 
tidal tail of the Arches cluster might be included in our 
outer cluster field. We take an average of control fields 
A, B, and C as the LF of the background population of 
the cluster. Note that this will probably overestimate the 
background population toward the Arches cluster slightly 
since we do not have a control field that is located at the 
opposite side of the control field B, and thus our estimate 
for the cluster LF/MF would be a lower limit. 

Figure [4] compares the dereddened ifg-band LFs for 
the background population (an average of control fields 
A, B, and C) and that for the intermediate region of 
the cluster field obtained with the medium field camera. 
The figure shows both raw and completeness-corrected 
LFs, and the latter of the cluster field is clearly larger 
than that of the background population at least down to 
Ks — 17 mag, which is the 50 % completeness limit for 
the cluster field. The 7?-band LFs give a similar result 
except that their 50 % completeness limit is 0.5-1 mag 
brighter than for the Kg band (the limit in H itself is 
1 mag brighter than that in Ks, but the typical apparent 
H — Kg colors of faint stars in the cluster field are 1.5-2 
mag). 

3.3. The Mass Function 

We convert the apparent magnitudes into masses us- 
ing the Geneva models. Here, we only use Ks-hand 
data as they provide deeper photometry. Figure [5] 
shows our completeness-corrected MF for the cluster in- 
termediate region after background subtraction, along 
with the MF from Stolte et al. (2005) that is also 
completeness-corrected, but not background-subtracted; 
both are scaled to an area of 5"-9" annulus. The MF 
of Stolte et al. (2005) shows a global turnover below 
logM/Mg ~ 0.35 (Af ~ 2.3 M©), which is also their 
50 % completeness limit, but our cluster MF, whose 50 % 
completeness limit is at logM/Mg = 0.1 (Af ~ 1.3 Mg), 
still has a significant amount of stars below M = 2.3 M© 
and keeps increasing at least down to M = 1 M© . Note 
that our MF increases even below our 50 % completeness 
limit. 

Stolte et al. (2005) claim that there may be a turnover 
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Fig. 5. — Background-subtracted mass functions derived from 
our ii's-band luminosity function for the 5"-9" annulus of the 
cluster field (solid) and the mass function of Stolte et al. (2005) 
for the 5"-10" annulus of the cluster (dotted). Error bars and the 
best-fit power-law relation (T = -0.91) for logM/M© = 0.1-1.7 
are shown for the former. The average of control fields A, B, and C 
were used as the background for the former, while the background 
population was not subtracted for the latter. Also plotted are the 
MFs for the same projected annulus from our Fokker-Planck calcu- 
lations that fit the observed MF at 2 Myr the best (offset by —0.5 
dex for clear presentation). The asterisks are for the calculation 
with nil = 0. 1 M0 and the diamonds are for m; = 1 Mq . Both cal- 
culations have an initial F of —1.1, and total mass of 4 X 10** Mq. 
The arrows indicate the locations of the 50 % completeness limit 
and the 6.3 Mq bump. 

near logM/Mg ~ 0.8 (M ~ 6.3 Mq) in their MF 
and that this might indicate a global decrease of the 
background-subtracted MF below that mass. However, 
our data, which have a quarter dex lower completeness 
limit in mass, indicate that the MF globally increases 
down to our completeness limit even after background 
subtraction. Therefore, we presume that the turnover 
around M = 6.3 Mq claimed by Stolte et al. (2005) is in 
fact a local bump in the MF as see in Figure [51 Eisen- 
hauer et al. (1998) find a similar bump in their LF of the 
Galactic starburst template NGC 3603, and show that 
such a bump is an indication of the young age (< 3 Myr) 
using the evolutionary pre-main-sequence tracks by Palla 
& Stabler (1993). 

When fit to a single power-law relation, our cluster 
MF gives a power-law exponent of F = —0.91 ± 0.08 
(the Sapleter slope is F = —1.35) for the mass range of 
logM/M© = 0.1-1.7. This is slightly steeper than the 
exponent for massive stars only (logM/MQ = 0.7-1.7), 
F = -0.71 ±0.15. 

^ These calculations consider internal dynamics of the cluster 
such as mass segregation and tidal evaporation among others; see 
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We have performed several Fokker-Planck calculations 
(for mi — 0.1 & IMq) and A^-body simulations (for 
mi — IMq) targeted for the Arches cluster with initial 
cluster conditions similar to those used in Kim, Morris, & 
Lee (1999) and Kim et al. (2000), one of which is no ini- 
tial mass segretation. We find that the F values for the 
projected annulus of 5"-9" increase by 0.1-0.2 during 
the lifetime of the Arches cluster, 2-2.5 Myr.^ There- 
fore, our results suggest that the IMF of the cluster has 
F = —1.0 ~ —1.1. Figure O shows two MFs for our annu- 
lus from our Fokker-Planck calculations that best match 
the observed MF at 2 Myr. These two calculations have 
considerably different mi's (0.1 & IMq), but result in 
nearly identical MFs for the annulus. Thus our estimate 
for the IMF from the simulations does not sensitively 
depend on the choice of mi . 

Using the red clump stars in the region \l\ < 2° and 
0°.5 < |6| < 1° that were observed with the Infrared 
Survey Facility telescope, Nishiyama et al. (2006) newly 
estimated the near-infrared extinction law toward the 
inner Galactic bulge. When assuming this extinction 
law, which is derived by averaging for rather larger re- 
gion of the GC, is applicable to the stars in our images, 
we find that the 50 % completness limit moves down to 
1 Mq, but there is almost no change in F for our annulus 
(-0.90 ±0.09). 
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